arXiv:1503.00236vl [quant-ph] 1 Mar 2015 


Enhancing parameter precision of optimal quantum estimation by direct quantum feedback 


Qiang Zheng lj 2 , Li Ge 2 , Yao Yao 2 , and Qi-jun Zhi 3 
1 School of Mathematics and Computer Science, Guizhou Normal University, Guiyang, 550001, China 
2 Laboratory for Quantum Optics and Quantum Information, 

Beijing Computational Science Research Center, Beijing 100084, China 
3 School of Physics and Electronics, Guizhou Normal University, Guiyang, 550001, China 

Various schemes have been proposed to overcome the drawback of the decoherence on quantum-enhanced 
parameter estimation. Here we suggest an alternative method, quantum feedback, to enhance the parameter 
precision of optimal quantum estimation of a dissipative qubit by investigating its dynamics of quantum Fisher 
information. We find that compared with the case without feedback, the quantum Fisher information of the 
dissipative qubit in the case of feedback has a large maximum value in time evolution and a smaller decay rate 
in the long time. 
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I. INTRODUCTION 

Quantum metrology H , which investigates the limit of the 
precision of parameter estimation bounded by the quantum 
mechanics, excites wide interest in recent years. One of the 
basic ideas is to surpass the shot-noise limit by making use 
of non-classical resources which provide, e.g., NOON state 
or entangled state. Quantum metrology has wide application 
in improving time and frequency standardsdetecting 
gravitational waves H , and magnetometry [6j], and so on. 

Quantum Fisher information (QFI) plays a central role in 
quantum metrology and quantum estimation theory. Accord¬ 
ing to quantum Cramer-Rao inequality, the QFI imposes an 
upper bound on the precision of parameter estimation. A 
larger QFI implies that the parameter can be estimated with 
a higher precision. In addition to the basic properties such 
as convexity and monotonicity, the QFI is also closely con¬ 
nected with other quantities, especially entanglement j7[], non- 
Markovianity fU], and spin squeezing [JSj, ITojl . More impor¬ 
tantly, the QFI closely connects with other branches of quan¬ 
tum physics, such as quantum phase transition ifllll . quantum 
clone [12], and quantum Zeno dynamics 0. 

Each realistic quantum system inevitably interacts with its 
environment, which leads to decoherence and considerable 
decrease of precision of parameter estimation 00. In or¬ 
der to enhance QFI, various methods, such as quantum er¬ 
ror correction, dynamical decoupling, decoherence-free sub¬ 
space, reservoir engineering, have been proposed to sup¬ 
press the decoherence. For instance, an approach based on 
decoherence-free subspace has been proposed j20ll to enhance 
the precision of estimation of atomic transition frequencies 
with ions stored in Paul traps subject to collective dephasing. 
Adopting the method of dynamical decoupling, the maximum 
Fisher information of a single qubit in a noisy system has been 
retrieved (21]. Inspired by Ref. [21], Tan et al. proposed a 
scheme to enhance the QFI of a qubit _by employing dynami¬ 
cal decoupling pulses (221. Berrada [|2J] and Chin et al. [ 24] 
also considered the QFI of a qubit subject to non-Markovian 
environment, which can be viewed as an engineered reservoir. 

Besides the above mentioned strategies, it is also interest¬ 
ing to use alternative methods to enhance QFI. Here we pro¬ 


pose using the quantum feedback to enhance the precision of 
parameter estimation of a dissipative qubit. Quantum feed¬ 
back mM has been considered as one of the essential ways 
to suppress the decoherence, which manipulates the system 
based on the information acquired by measurement. Com¬ 
pared with the other control approaches, it operates simply by 
feeding back some information, which is proportional to the 
measurement signal, directly to the controlled system [27]. 
Generally, it is well known that the feedback can only partly 
compensate the decoherence effect of an open system gov¬ 
erned by the Lindblad master equation. In the method of di¬ 
rect feedback, there are various schemes to measure the out¬ 
put field, such as homodyne or direct photodetection. We will 
adopt the latter scheme and only consider the perfect detection 
efficiency in this paper. 

As a first step to explore the possibility that quantum feed¬ 
back can beat the shot-noise limit and eventually achieve the 
Heisenberg limit, we study a qubit system. Moreover, moti¬ 
vated by recent experimental achievement on quantum feed¬ 
back control on a qubit system f28l l29h . in this paper we 
study the parameter estimation of a dissipative qubit in a cav¬ 
ity which provides the quantum feedback to the qubit. For 
the purpose of exploring the effect of quantum feedback, we 
only focus on the QFI of the qubit with considering the opti¬ 
mal conditions for saturated quantum Cramer-Rao inequality. 
Under the effect of quantum feedback, we find that: (i) in the 
long-time limit, the decay rate of the QFI with feedback is 
slowed down; and (ii) the maximum value of the QFI in the 
evolution is also enhanced. These results imply that the quan¬ 
tum feedback can enhance the precision of parameter estima¬ 
tion. In addition, we find there is an optimal external driving 
strength for the qubit maximizing the QFI of the steady state. 

This paper is organized as follows. The basic properties of 
the QFI are reviewed in Sec. II. In Sec. Ill, we investigate the 
feedback effect on the QFI of a qubit interacting with a cavity, 
and find that the QFI is improved by feedback and discuss 
the behaviours of the QFI in the steady state with the external 
driving of the qubit. Finally, a summary is provided in the last 
section. 
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II. QUANTUM FISHER INFORMATION 


III. ENHANCING QFI OF A QUBIT WITH FEEDBACK 


In this section, we will review the main aspects of QFI. Let 
ip denote a single parameter to be estimated, and Pi(ip) be 
the probability density with measurement outcome {xi} for a 
discrete observable X conditioned on the fixed parameter ip. 
The classical Cramer-Rao inequality 13011 gives the bound of 
the variance Var(</3) for an unbiased estimator p> 

Var (<p) > ( 1 ) 

where the classical Fisher information is defined as Mh v = 

Extending to quantum regime, in order to determine the ul¬ 
timate bound to precision posed by quantum mechanics, the 
Fisher information must be maximized over all possible mea¬ 
surements Q. By introducing the symmetric logarithmic 
derivative L v determined by = ^{p v L v + L v p ip ), the 
so-called quantum Cramer-Rao inequality gives a bound to 
the variance of any unbiased estimator 0311 : 

Var (^) > _i_ > ^. ( 2 ) 

It is worth mentioning that the symmetric logarithmic deriva¬ 
tive L v provides an optimal measurement, that is, the inequal¬ 
ity in Eq. (0 can be saturated via using a measurement com¬ 
posed by the eigenvectors of L v . Here, the QFI of a quantum 
state p v with respect to the parameter ip is defined as 0411 

F v = Tr { P(p L 2 v ). (3) 

With the spectrum decomposition p^, = Afc|fc)(fc|, the 
QFI is completely determined by 


A. Model 

We consider a system consisting of a qubit resonantly in¬ 
teracting with a single-mode cavity, as shown in Fig. [T| The 
(real) coupling strength between the qubit and the cavity is g, 
the cavity mode a is damped at rate n, the two levels of the 
qubit are |e) and |p), and its spontaneous emission rate of the 
excited state is 70 . Without feedback, the time evolution of 
the whole system is described by the Lindblad master equa¬ 
tion (h = l) 

+ KD[a\p + ^oD[a-]p, ( 8 ) 

where the superoperator is defined as V[c]p = cpc* — A(ct cp+ 
pc^c), and the Hamiltonian of the model under consideration 
is given as 

H = §07 + g(a + a + h.c.). (9) 

Here we have assumed that the qubit is subject to an external 
classical drive with strength O. The Pauli operators a z = 
|e)(e| - \g)(g\ and a x = ct_ + a + with <7_ = |p)(e| and 
<r+ = |e)(p| being the lowering and raising operators of the 
qubit, respectively. 

Under the condition that the cavity decay k is much larger 
than the other relevant frequencies of the system, the cavity 
mode can be adiabatically eliminated and the effective damp¬ 
ing rate of the qubit is obtained as 0 

7 = g 2 / K - (io) 


K 


E 

fc,Afc>0 


(dyAfc ) 2 
A k 


k , 


E 

k r ,Afc -t-A^./ >0 


2(Afc — Afc') 2 
Afc + Afc' 


\(k\d ip k')\ 2 . 


(4) 


The first term in this equation is just the classical Fisher infor¬ 
mation, and the second term can be considered as the quantum 
contribution. And the symmetry logarithmic derivative can be 
obtained as 


L(p — 


E 


k , k r , Afc+A fc />0 


2(k\d v p\k') 

Afc + Afc' 


1 km. 


(5) 


Usually it is difficult to give the QFI in the form the ma¬ 
trix density for a general system. Fortunately, the QFI of 
the two-dimensional density matrix has obtained explicitly in 
Refs. JMEU] as 

F<p = Tr [{d v p) 2 ) + d^Ti -[{pd^p) 2 }. ( 6 ) 

Moreover, the QFI is also related to the Bures distance lf33ll 
through 

[Pv > P<p+d(p] = ~^F'<pdtp~, (7) 

where the Bures distance which measures the distance be¬ 
tween two quantum states p and a is defined as 0 
DbIp,*} = [2(1 - TrvV/Vpi/ 2 )] 1 / 2 . 


In what follows, we will consider the limit 7 70 such that 

the spontaneous emission of the qubit can be omitted. 

Now let us consider the system subject to the feedback, as 
shown in Fig. [Q The output from the cavity is measured by a 
detector D, then the signal I(t) from the detector D triggers 
the control Hamiltonian iFfb = I(t)B. For a feedback scheme 
based on photodetection measurement, the unconditional state 
of the qubit is described by li38il39ll 

^ =-iln[a x ,p\+'yV[Ua-\p. (11) 

Rewriting the last term of Eq. (fTTb explicitly: V[Ua-]p = 
Ucj-pa + U * — A(cr + < 7 _p + pcr + (j_), the effect of the feed¬ 
back is obvious. Conditioned on a detection event cr_pcr + , 
the unitary evolution U = e~ lB6t is immediately applied on 
the system lasting finite time St. Throughout this paper, the 
other relevant parameters are scaled by St, such as that t or 
7 denotes the scaled time T or decay rate 7 St, and we set 
St = 1. Taking the constraint UW = 1 into consideration 
(which is required by experiments), we choose the feedback 
operator as 14011 

U = e= cos(A)/ 2 + i sin( A) [a x sin(/3) + a y cos(/3)]. 

112 ) 

Here I 2 is 2 x 2 identity operator, = (<r x ,a y ,a z ) and A = 
(Asin(/3),Acos(/3),0). The two parameters A and /? fully 
characterize the feedback operator U. 
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FIG. 1: (Color online) Schematic diagram of the parameter estima¬ 
tion with feedback. The qubit is coupled to the cavity resonantly. 
The output of the cavity is measured, and the measurement is used to 
trigger a feedback Hamiltonian. 


B. Without feedback 


In this subsection, we investigate the dynamics of the QFI 
in absence of the feedback. Here we focus on that the external 
driving is switched off, the case fl 7 ^ 0 will be discussed in 
Sec. Ill D. The effective damping in Eq. (ITT1) can be consid¬ 
ered as a quantum channel for the qubit. Then we adopt the 
quantum parameter estimation theory to estimate the effective 
damping parameter 7 . 

For an superposition initial state |t/>o) = (cos($)|e) + 
sin(r?)|g))/, we first concentrate on 1 } = 7 t/ 4 . The effect 
of population on |e) will be discussed later. For this equal- 
weighted superposition initial state, the evolved density ma¬ 
trix of the qubit can be exactly solved, which is given as 

*>-(aw wife)) 

with the elements 

Pit, = e _ * 7 /2, p 12> off (t) = e- t 7 / 2 /2, (14) 


under the condition without the feedback, i.e., by setting A = 

7T. 

Adopting Eq. (O by some straightforward calculations, the 
QFI of p{t) in Eq. ( IT3l ) with respect to 7 is obtained as 


F^t) 


i 2 e-‘1'(2e tl '-l) 

4(e t T-l) 


(15) 


In the limit t —> 0, we find that 



(16) 



t 



FIG. 2: (Color online) The evolution of the eigenvalue of the qubit 
density matrix under the effect of the control parameters A in panel 
(a) and /3 in panel (b). Here /3 = 0 in (a) and A = 7t/3 in (b), re¬ 
spectively. The other parameter is 7 = 0.1, and A = ix corresponds 
to turning off the quantum feedback. All the relevant quantity, scaled 
by the feedback lasting time St, is dimensionless mentioned here and 
in the other figures (setting St = 1). 


C. Enhancing QFI by feedback 

In order to show the effect of quantum feedback, in Fig. [2] 
we plot the evolution of the smaller eigenvalue of the qubit 
density matrix p(t) under the effect of the feedback parame¬ 
ters. Fig. [3a) shows that the (smaller) eigenvalue of the qubit 
is increased by the feedback and has the slow decay in the 
long-time limit with A = 7t/3. Fig. |2jb) presents that the 
eigenvalue decreases with the increase of /3 in the initial time 
and is almost independent of the feedback parameter 8 in the 
long-time limit. Thus, we keep (3 = 0 in the following discus¬ 
sions as a good approximation. 

When the feedback is turned on (but still in absence of qubit 
driving), the evolved elements of the qubit density matrix are 
given as 


This implies that the QFI increases with time in the initial 
evolution. On the other hand, in the long time limit t —> 00 , 
the QFI reduces to 

F 7 (t) f 2 e _ 7 F */2, (17) 

which exponentially decays to zero with rate jp = 7 . The 
similar result has also been obtained in Ref. 14 ill , although 
a different estimated parameter was considered in Ref. [41], 
i.e., a phase parameter within input states, which experience 
amplitude damping or spontaneous emission before measure¬ 
ment. 


Pu(t) — e tlq /2 , pi 2 (i) — P 12 , o//(l + y), (18) 

resulting from Eq. (fTTl) with 8 = 0. Here, it is apparent that 

7 , = 7cos(A) 2 , Y = (1 — e -t7COS ( 2j4 )/ 2 ) tan(2A), 

(19) 

are attenuating factors of qubit density matrix induced by the 
quantum feedback. 

For general feedback parameter A, the analytical expres¬ 
sion of the QFI is very cumbersome. As a warm-up, we study 
two relatively simple situations. For A = 7 r/ 2 , the effect of 
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FIG. 3: (Color online) The evolution of the QFI F 7 (/,) with different 
decay rates 7 = 0.1 (a) and 7 = 0.5 (b). The dashed line (A = 7t/3) 
correspond to turning on the feedback, while the solid line (A = n) 
to off the feedback, respectively. 


the feedback is equivalent to qubit dephasing: only the off- 
diagonal elements decays exponentially with a rate 7 q =7/2. 
In this case, the QFI is obtained as 

Fy(t) = 4 ( e t 7 _ 1 ), (20) 

which also decays exponentially to zero with rate 7 ^ = 7 in 
the long-time limit. Therefore, there should exist an optimal 
time to maximize the value of the QFI. For A = 7 r/ 3 , which 
makes all the elements of the qubit density matrix in Eq. (fTSl i 
exponentially decay to zero with rates = 7 / 4 . In the long¬ 
time limit, the QFI is given as 

F 7 (t) ~ ^e-^t 2 . ( 21 ) 

with a rate 7 ^ = 7 / 4 . This smaller decay rate means the pre¬ 
cision of parameter estimation is enhanced by the feedback. 
In Fig. [3] we plot the QFI with respect to time t. This figure 
shows that in the short time, the value of the QFI without feed¬ 
back is larger than that with feedback. Fig.[3]also displays that 
the QFI with the feedback decays slowly than that without the 
feedback in the long time. Moreover, the maximum value of 
the QFI is enhanced by the feedback. With the increase of the 
decay rate, the value of the QFI drops considerably. 

For the general parameter, the QFI at the fixed time is plot¬ 
ted in Fig. [4] in order to give further information of the feed¬ 
back effect. This figure is consistent with Fig. [3] in the long¬ 
time limit, the QFI is enhanced by the feedback. The peak 
structure in Fig.@]implies that the estimation precision can be 
strengthened by carefully tuning the feedback strength A to an 



FIG. 4: (Color online) The fixed time QFI F 7 (fo) in terms of the 
feedback parameters A with /3 = 0. The dashed (solid) line corre¬ 
sponds to the time to = 20 (to = 300), respectively. The feedback 
is switched off at A = 0, 7 r. The other parameters are 7 = 0.1. 


40 r 
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FIG. 5: (Color online) The variation of of maximum value of the 
QFI Fm with respect to A. It is periodic function of A with a period 
A = 7 r. The parameter is chosen as 7 = 0.1. 


optimal value. We also find that the QFI is a periodic function 
of A with period T = tt. 

As the maximum value of QFI implies the largest precision 
to estimate a parameter, in Fig.[5]we plot Fm as a function of 
A, with Fm defined as 


F m = max[F 7 (£)]. ( 22 ) 


Without feedback (A = 0, 7 r), the maximum value of the QFI 
Fm ~ 29.0. At presence of the feedback, Fm is enhanced, 
especially in the neighborhood of A = tz/?> , Fm ~ 38. This 
figure also shows that near A = 7 t/ 2 , the quantum feedback 
paly a negative role for the QFI. 

We have investigated the QFI for a specific initial state: an 
equal-weight superposition pure state. Next, we will study the 
behaviors of the QFI for other initial states, especially a mixed 
one. For an initial mixed state po = er|e)(e| + (1 — e)\g)(g\ 
(e £ [0,1)), the QFI for the general feedback parameters A 
and /? is also very complicated. For the feedback value A = 
7 t/3 , the QFI is obtained as 


Fy(t) 


m(t)z + V2(t)£ 2 + 7) 3 ft)£ 3 

m(t)£ + r]5 


(23) 
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with 

771(4) = —e 3 * 7,/2 4 2 /2, 772(4) = (—12c 3 * 7 / 4 + 6e* 7 )4 2 , 

773 ( 4 ) = 6e* 7 ,/ 2 4 2 , 775 ( 4 ) = — 8e 7 * 7 / 4 

774 ( 4 ) = 8(e 3 * 7 / 2 + 3(e‘ 7 - 2e 5 * 7 / 4 + e 3 ‘ 7 /2)). 

(24) 

In the limit 4—^0, this QFI reduces to 

i^(4) =C(e)4 2 (25) 

with an enhancing factor £(e) = being a 

monotonically increasing function of e. This means that the 
large probability in the excited state can enhance the precision 
of parameter estimation. In the long-time limit, 

F 7 (4) « jg£e -t7 / 4 4 2 , (26) 

which also goes up with the increase of e. 

The effect of the initial state on the QFI is further studied 
by the numerical simulation. We find that the maximum value 
of the QFI linearly (approximately) increases with the prob¬ 
ability in the excited state s. The physical reason of these 
results can be understood as follows: The information of 7 
comes from the effective damping of the excited state. The 
initial state with a high probability in the excited state should 
produce a larger QFI value. In absence of the feedback (i.e., 
A = 7 r), the QFI is gained as Fly (4) = e t*_ e , which has larger 
exponential decay than that of A = 7t/3 in the long-time limit. 
As a short summary, the QFI can be enhanced by the feedback, 
independent of the input state. 


D. The QFI of the steady state with driving 

In the previous subsection, we study the QFI with respect 
to the effective damping parameter 7 with switching off the 
external driving. The main result is the QFI can be enhanced 
by the feedback. However, in this case of no external driving, 
the QFI decays to zero exponentially in the long-time limit. 
Here we extend our study to the case of the external driving 
strength 11/0. In this situation, the time-evolving density 
matrix of the qubit at any time can not be analytically given. 
However, one can obtain the non-zero steady state in the long¬ 
time limit (actually, in many realistic information processing, 
the steady state is more useful) as 

Ps, 11 = fl 2 /M, p s , 12 = N/M (27) 

with 

N = Q[f2sin(2A) cos(/3) — icos(A) 2 7 ], 

(28) 

M = Q, sin(2A)7 sin(/3) + 2 Q 2 + cos(A) 2 7 2 . 

We have checked that the QFIM of this steady state satisfies 
Det(F(< 7 , k)) = 0. The above expressions for the steady state 
is also complicated, then we investigate its property by con¬ 
sidering specific parameters. In the strong-driving limit, i.e., 
fl is much larger than the other parameters, the elements of 



Q. 

FIG. 6 : (Color online) The QFI of the steady state with respect to 
f i. There is an optimal driving strength to maximize the F 7 (og). 
The other parameters are A = n/3. All the relevant quantity is 
dimensionless, such as H is scaled by SI St. 


the state reduce to p s ,n = 1/2 and p s ,i2 = 0 , and the corre¬ 
sponding QFI with respect to 7 equals to zero. 

By setting A = n/3, the elements of the density matrix 
translate to p s ,n = Q 2 /x, p s ,i2 = fi(\/3^ — *7/2)/%, with 
X = 7 ) 7 2 + 4fl 2 . The QFI of this steady state with respect to 
7 is obtained as 


? _ i6n 2 (37 2 +n 2 ) 

7 — (3^q^n4jFP+8TFF' 


(29) 


It is obvious that without the driving (Q = 0), the QFI of 
the steady state equals to zero. In the limit fl —> 0, the QFI 
becomes F 7 oc SI 2 , which goes up with the increase of Q. 
Taking into consideration the fact that the QFI is close to zero 
in the large driving limit, there should have an optimal driving 
strength (),,, to maximize the QFI. In Fig. [ 6 ] we show the vari¬ 
ation of the QFI with respect to f2, which is consistent with 
the analytical analysis. This figure also shows that with the 
increase of 7 , the value of the QFI drops considerably, which 
implies the precision of estimation of 7 decreases. 


IV. CONCLUSION 

In summary, in the system that the qubit couples to the 
cavity, we demonstrate that the quantum feedback can 
enhance the the precision of parameter estimation. Quantum 
Fisher information is used as a measure to describe the 
parameter-estimation precision as one can use the observable 
composed by the eigenvectors of the symmetry logarithmic 
derivative to saturate the quantum Cramer-Rao inequality. 
Different from the other quantum correlation such as quan¬ 
tum entanglement and quantum steering, the quantum Fisher 
information characterizes the sensitivity of a quantum state 
with respect to the change of a certain parameter. We find that 
in presence of the feedback, the decay rate of the quantum 
Fisher information is slowed down and the maximum value 
of the quantum Fisher information is also strengthened. In 
future, it’s an interesting topic to explore the effects of a 
homodyne measurement scheme and the detection efficiency 
on the quantum Fisher information. In addition, the quantum 
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Fisher information of multi-qubits PkZjl are another valuable 
subject. 
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